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Abstract. A detailed study of the modified gravity, f(R) models is performed, using that the 
Noether point symmetries of these models are geometric symmetries of the mini superspace 
of the theory. It is shown that the requirement that the field equations admit Noether point 
O ■ symmetries selects definite models in a self-consistent way. As an application in Cosmology we 

consider the Friedman -Robertson- Walker spacetime and show that the only cosmological model 
, which is integrable via Noether point symmetries is the (i? b — 2A) model, which generalizes 

the Lambda Cosmology. Furthermore using the corresponding Noether integrals we compute 
the analytic form of the main cosmological functions. 
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1. Introduction 

The recent cosmological data indicate that the universe is spatially flat and has suffered two 
acceleration phases. An early acceleration phase (inflation), which occurred prior to the radiation 
CN ■ dominated era and a recently initiated accelerated expansion. 

An easy way to explain this expansion is to consider an additional fluid with negative equation 
of state parameter, usually called dark energy, that dominates the universe at late times. In 
spite of that, the absence of a fundamental physical theory, regarding the mechanism inducing 
the cosmic acceleration, has given rise to a plethora of alternative cosmological scenarios. Most 
of them are based either on the existence of new fields in nature (dark energy) or in some 
modification of Einstein's general relativity (GR), with the present accelerating stage appearing 
as a sort of geometric effect (" geometrical" dark energy) . 

The simplest dark energy probe is the cosmological constant A (vacuum) leading to the 
ACDM cosmology [TH3]. However, it has been shown that ACDM cosmology suffers from two 
major drawbacks known as the fine tuning problem and the coincidence problem [5J. Besides 
ACDM cosmology, many other candidates have been proposed in the literature, such as time- 
varying A(i) cosmologies, quintessence, k— essence, tachyons, modifications of gravity, Chaplygin 
gas and others [5HT0]. 

There are other possibilities to explain the present accelerating stage. For instance, one may 
consider that the dynamical effects attributed to dark energy can be resembled by the effects of a 
nonstandard gravity theory. In other words, the present accelerating stage of the universe can be 
driven only by cold dark matter, under a modification of the nature of gravity. Such a reduction 
of the so-called dark sector is naturally obtained in the f(R) gravity theories In the original 
nonstandard gravity models, one modifies the Einstein-Hilbert action with a general function 



f(R) of the Ricci scalar R. The f(R) approach is a relative simple but still a fundamental tool 
used to explain the accelerated expansion of the universe. A pioneering fundamental approach 
was proposed long ago with f(R) = R + mR 2 [12]. Later on, the f(R) models were further 
explored from different points of view in [13H15] and indeed a large number of functional forms 
of f(R) gravity is currently available in the literature |16H21j . 

The aim of the present work is to investigate which f(R) models admit extra Noether point 
symmetries and use the first integrals of these models to determine analytic solutions of their 
field equations. The idea to use Noether symmetries in cosmological studies is not new and 
indeed a lot of attention has been paid in the literature (see |22H37j ). 

The main reasons for the consideration of this hypothesis is that (a) the Noether point 
symmetries provide integrals, which assist the integrability of the system, (b) is a geometric 
criterion because the Noether symmetries associated with the geometry of the field equations. 
A fundamental approach to derive the Noether point symmetries of a given dynamical system 
moving in a Riemannian space has been proposed in [38]. A similar analysis can be found 
in [31110]. 

The structure of the paper is as follows. The basic theoretical elements of the problem are 
presented in section [21 where we also introduce the basic FRW cosmological equations in the 
framework of f(R) models. The geometrical Noether point symmetries and their connections 
to the f(R) models are discussed in sections [5j In section [6] we provide analytical solutions for 
those f{R) models which are Liouville integrable via Noether point symmetries. In section [7] 
we study the Noether symmetries in spatially non flat f(R) cosmological models. Finally, we 
draw our main conclusions in section 



2. Cosmology with a modified gravity 

Consider the modified Einstein-Hilbert action: 



S 



d x^J—g 



(1) 



where C m is the Lagrangian of dust-like {jp m 
with respect to the metric^] we arrive at 



= 0) matter and k 2 = 8irG. Varying the action 
2Df - (/ - Rf 



6? t = k 2 TH 



(2) 



where the prime denotes derivative with respect to R, Gy is the Einstein tensor and Ty is the 
ordinary energy-momentum tensor of matter. Based on the matter era we treat the expanding 
universe as a perfect fluid which includes only cold dark matter with comoving observers 
XJ^ = Sq. Thus the energy momentum tensor becomes T^ v = p m Uf_iU u , where p m is the energy 
density of the cosmic fluid. 

Now, in the context of a flat FRW model the metric is 



ds 1 



-dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). 



(3) 



The components of the Einstein tensor are computed to be: 

G° = -3H 2 , G a b = -<5£ (2H + 3# 2 ) . (4) 



1 We use the metric i.e. the Hilbert variational approach. 



Inserting into the modified Einstein's field equations ([2]), for comoving observers, we derive 
the modified Friedman's equation 



3/'tf 2 = k 2 p m + f -^-l - 3Hf"R (5) 



2f'H + 3f'H* = -2Hf"R - ( /"' if + f" R ) - f -J^L. ( 6 ) 



The contraction of the Ricci tensor provides the Ricci scalar 

R = sTRp, = 6 + ^) = 6(2^ 2 + H) . (7) 

The Bianchi identity T^ u = leads to the matter conservation law: 

Pm + 3Hp m = (8) 

whose solution is 

Pm = Pm,oa~ 3 . (9) 

Note that the over-dot denotes derivative with respect to the cosmic time t and H = a/a is the 
Hubble parameter. 

If we consider f(R) = R then the field equations (|2J) boil down to the Einstein's equations 
a solution of which is the Einstein de Sitter model. On the other hand, the concordance A 
cosmology is fully recovered for f(R) = R — 2 A. 

From the current analysis it becomes clear that unlike the standard Friedman equations in 
Einstein's GR the modified equations of motion (J5J) and © are complicated and thus it is 
difficult to solve them analytically. 

We would like to stress here that within the context of the metric formalism the above f(R) 
cosmological models must obey simultaneously some strong conditions [42J. These are: (i) / > 
for R > i?o > 0, where Rq is the Ricci scalar at the present time. If the final attractor is a 
de Sitter point we need to have / > for R > Ri > 0, where Ri is the Ricci scalar at the 
de Sitter point, (ii) /" > for R > R > 0, (iii) f(R) ra R - 2A for R > R and finally (iv) 

< H-(r) < 1 at r = -Sj- = -2 



3. Modified gravity versus symmetries 

In the last decade a large number of experiments have been proposed in order to constrain 
dark energy and study its evolution. Naturally, in order to establish the evolution of the dark 
energy ("geometrical" in the current work) equation of state parameter a realistic form of H(a) 
is required while the included free parameters must be constrained through a combination of 
independent DE probes (for example SNIa, BAOs, CMB etc). However, a weak point here is the 
fact that the majority of the f(R) models appeared in the literature are plagued with no clear 
physical basis and/or many free parameters. Due to the large number of free parameters many 
such models could fit the data. The proposed additional criterion of Noether point symmetry 
requirement is a physically meaning- full geometric ansatz. 

According to the theory of general relativity, the space-time symmetries (Killing and 
homothetic vectors) via the Einstein's field equations, are also symmetries of the energy 
momentum tensor. Due to the fact that the f(R) models provide a natural generalization 
of GR one would expect that the theories of modified gravity must inherit the symmetries of 
the space-time as the usual gravity (GR) does. 



Furthermore, besides the geometric symmetries we have to consider the dynamical 
symmetries, which are the symmetries of the field equations (Lie symmetries). If the field 
equations are derived from a Lagrangian then there is a special class of Lie symmetries, the 
Noether symmetries, which lead to conserved currents or, equivalently, to first integrals of the 
equations of motion. The Noether integrals are used to reduce the order of the field equations or 
even to solve them. Therefore a sound requirement, which is possible to be made in Lagrangian 
theories, is that they admit extra Noether symmetries. This assumption is model independent, 
because it is imposed after the field equations have been derived, therefore it does not lead to 
conflict with the geometric symmetries while, at the same time, serves the original purpose of 
a selection rule. Of course, it is possible that a different method could be assumed and select 
another subset of viable models. However, symmetry has always played a dominant role in 
Physics and this gives an aesthetic and a physical priority to our proposal. 

In the Lagrangian context, the main field equations ([5]) and ©, described in section [2j can 
be produced by the following Lagrangian: 

L (a, a, R, r) = 6a/' d 2 + 6a 2 /" hR + a 3 (f'R - f) (10) 



in the space of the variables {a, R}. Using eq. (|10p we obtain the Hamiltonian of the current 
dynamical system 

E = 6a/' d 2 + 6a 2 /" aR — a 3 (^f R — A (11) 

or 

* ' (12) 



E = 6a 3 



f'H 2 --^ 7 [{fR-f)-QRHf 



6/ 

Combining the first equation of motion ([5]) with eq. (|12p we find 



The latter equation together with p m = p m §a 3 implies that 

E E 

PmO = ^Tf ^ ^mPcrfi = WT^ => E = QVl m H Q (14) 

where O m = p m o/p C r,o, Pcr,o = 3Hq/Ic 2 is the critical density at the present time and Hq is the 
Hubble constant. 

We note that the current Lagrangian eq. (|1U|) is time independent implying that the dynamical 
system is autonomous hence the Hamiltonian E is conserved (^ = 0). 



4. Noether symmetries 

Before we proceed we review briefly the basic definitions concerning Lie and Noether point 
symmetries of systems of second order ordinary differential equations (ODEs) 

■£ = J (t,x j ,x j ) . (15) 

The one point parameter transformation 



t = t + e£(t,j) (16) 
= x i + erf{t,x i ) (17) 



with generator X = £ (t, a^) dt + f? 4 (t, x) di is a Lie point symmetry of the system of ODEs (|15p 
if the following condition is satisfied [39^40] 



xW(x i -u(t,x j ,x j ))=0 (18) 
where X^ is the second prolongation of X defined by the formula 

X® = + 17*$ + (if - x't) d^i + (ff - tf'i - 2x i i) d&. (19) 



Condition ([18]) is equivalent to the relation 

XW,A] = A(x a ),4 (20) 

where XW is the first prolongation of X and A is the Hamiltonian vector field 

A = d t + xd x + J (t, x j , x j ) d ±l . (21) 

If the system of ODEs results from a first order Lagrangian L = L [t,xi,x 3 } , then a Lie 
symmetry X of the system (|15|) is a Noether symmetry of the Lagrangian if the additional 
condition is satisfied 

x[1]L dJ L = dg 

dt dt y ' 

where g = g (t, x 3 ) is the gauge function. To every Noether symmetry there corresponds a first 
integral (a Noether integral) of the system of equations (|15p which is given by the formula: 

BT 

I = ZE H -— i r, i + g (23) 
where Ejj is the Hamiltonian of the Lagrangian 

BT 

= (24) 
The vector field X in the augmented space {t, a, R} is 

X = £(t, a, R) 8 t + r?« (t, a, R) 8 a + t}<® (t, a, R) 8 R (25) 
and the first prolongation 

XW = £dt + rPd a + V {2) d R + (n® - di) da + (fi<?> - ItA Br. (26) 



Having given the basic formula for the Noether symmetries we look for analytic solutions of 
the dynamical system with Lagrangian (|1U|) with the use of Noether Integrals. 

5. Noether symmetries of / (R) gravity 

The Noether condition (|22p for the Lagrangian fllQ|) is equivalent with the following system of 
eight equations 



e,« = o 
i R = o 



(27) 
(28) 



a 2 f"v ( S = (29) 

f v m + afrjW + 2a /'r ? ( a 1 ) + a 2 f"r,^ - ~a/'0 = (30) 

2af" V W + a 2 r?? (2) + a 2 fv (l) + 2a/7? (l) + a 2 A g _ l a 2 f% = (31) 

- 3a 2 Rf'r,^ + 3a 2 /r/ (1) - a 3 Rf" V ^ + a 3 (/ - /i?) ^.t + <?,* = (32) 

Uaf'r,^ + 6a 2 f"^ + a 3 (/ R ~ f) £,a - 9,a = (33) 

GaY'rA 1 ) + a 3 f/'i? - f) in - g, R = (34) 



The solution of the system (|27p -(|34 p will determine the Noether symmetries. 

Since the Lagrangian (|10p is in the form L = T (a, a, R, Rj — V (a, R), the results of |38j can 

be usedE). The kinematic term defines a two dimensional metric in the space of {a, R} with line 
element 

ds 2 = Uaf'da 2 + 12a 2 f"da dR (35) 

while the "potential" is 

V(a,R) = -a 3 (f'R- f) . (36) 

The Ricci scalar of the two dimensional metric (|35p is computed to be R = 0, therefore the 
space is a flat spac^l with a maximum homothetic algebra. The homothetic algebra of the 
metric (1331) consists of the vectors 



f I 1 f 

K 1 = ad a - 3—8 R , K 2 = -d a - -2jjjd R 



a f" ' 2 2/ 

where K are Killing vectors (K 2,3 are gradients) and H is a gradient Homothetic vector. 
Therefore applying theorem 2 of [38] we have the following cases: 

Case 1: If / (i?) is arbitrary the dynamical system admits as Noether symmetry the X 1 = 
dt with Noether integral the Hamiltonian E. 

3 

Case 2: If / (R) = R2 the dynamical system admits the extra Noether symmetries 

X 2 = K 2 , X 3 = tK 2 (37) 

X 4 = 2td t + H+^K l . (38) 



with corresponding Noether Integrals 

dt 



I 2 = ^L[aVR) (39) 



J 3 = t j t ^aVR^j - aVR (40) 



2 Where T is the "kinetic" term and V is the "potential 
All two dimensional Riemannian spaces are Einstein sp 
symmetric [41] and if R — 0, the space admit gradient Homothetic vector, i.e. is flat. 



3 All two dimensional Riemannian spaces are Einstein spaces implying that if R = const the space is maximally 



_ fl 3 

h = 2tE - 6a 2 aVR - 6-=R. (41) 



'R 

the non vanishing commutators of the Noether algebra are 

[X\X 3 ]=X 2 [X\X A ]=2X 1 

[x 2 ,x*] = 8 -x 2 [x 3 ,x'] = 2 -x 3 

7 

Case 3: If / (R) = R» the dynamical system admits the extra Noether symmetries 

X 5 = 2td t + H , X 6 = t 2 d t + tli (42) 
with corresponding Noether Integrals 

'-**-!!(■*-') (43) 

h = t'E - (a?nri) + 1 a 3 nri. (44) 

and the non vanishing commutators of the Noether algebra are 

[X\X 5 ]=2X l [X 1 ,X 6 ]=X 5 [X 5 ,X 6 ]=2X 6 

From the time dependent integrals (|43p . (|44p and the Hamiltonian we construct the Ermakov- 
Lewis invariant 13.1 11] 

£ = 4/ 6 £-/ 2 (45) 

3 

Case 4: If / (R) = (R — 2A) 2 the dynamical system admits the extra Noether symmetries 

X 2 = e^K 2 , X 3 = e-^K 2 (46) 
with corresponding Noether Integrals 

h = e^ 1 (jt (ay/R - 2L) - 9y/may/R - 2A^j (47) 

J 3 = e -V^* M ^Vi? - 2l) + 9yftHay/R - 2A^j 
where m = |A. The non vanishing commutators of the Noether algebra are 



(48) 



[x\x 2 ] = ^x 2 [x 3 , X 1 ] = ^Rx 3 

From the time dependent integrals (|47p . (|48p we construct the time independent integral 
^23 = hh- 

7 

Case 5: If / (R) = (R — 2A) 8 the dynamical system admits the extra Noether symmetries 

X 5 = ^=e 2 ^ l d t + e 2 ^ lt H (49) 
X 6 = ^—e^^dt + e" 2 ^* H (50) 



with corresponding Noether Integrals 



2-E-~( a 3 (R - 2A)' 1 A + ^^a 3 (R ~ 2A)~B 



m 8 dt \ 7 4 



/m 8 dt \ y ' ) 4 

and the non vanishing commutators of the Noether algebra are 

[X\X 5 ] = 2 V / ^X 5 [X^X 1 ] = 2^X 



(51) 
(52) 



From the time dependent integrals (|43p . (|44p and the Hamiltonian we construct the Ermakov- 
Lewis invariant |44j 

4> = E 2 -hh (53) 

Case 6: If / (R) = R n (with n ^ 0, 1, |, |) the dynamical system admits the extra Noether 
symmetry 

I^^ + H+^-^K 1 (54) 
with corresponding Noether Integral 

I 7 = 2tE - 8na 2 R n ~ 1 a (2 - n) - 4na 3 R n ~ 2 R {2n - 1) (n - 1) . (55) 

and the commutator of the Noether algebra is [X , X 7 ] = 2X 1 . 

We note that the Noether subalgebra of case 2, {X 1 , X 2 , X 3 } and the algebra of case 4 
{AT 1 , X 2 , X 3 } is the same Lie algebra but not in the same representation. The same observation 
applies to the subalgebra of case 3 {X 1 , X 5 , X 6 } and the algebra of case 5 {X 1 , X 5 , X 6 }. This 
connection between the Lie groups is useful because it reveals common features in the dynamic 
systems, as is the common transformation to the normal coordinates of the systems. 

For the cosmological viability of the models see 



6. Analytic Solutions 

Using the Noether symmetries and the associated Noether integrals we solve analytically the 
differential eqs.([5]), Q and ([7]) for the cases where the dynamical system is Liouville integrable, 
that is for cases 2-5. Case 6 (i.e. / (R) = R n ) is not Liouville integrable via Noether point 
symmetries, since the Noether integral (|55|) is time dependent^. 

6.1. Power law model R 11 with fx = | 

In this case the Lagrangian eq. (|10p of the f(R) = R* model is written as 

Q^2 (Tj^ 3 

L = 9aV~Rd 2 + — =aR + — i?2 (56) 
2VR 2 

Changing the variables from (a, R) to (z, w) via the relations: 

"=(!)"' ^ fl = T (57) 

4 In the appendix |Appcndix A| we present special solutions for the / (R) = R n model, using the zero order 
invariants. 



the Lagrangian (j56j) and the Hamiltonian (jlip become 

L = zw + V w 3 (58) 

E = zw- V w 3 (59) 
where Vq = ^. The equations of motion in the new coordinate system are 

w = (60) 
z-3\W = (61) 

The Noether integrals (|39p . (|40p in the coordinate system {z,y} are 

I x = w , I 2 = tw — w (62) 

The general solution of the system is: 

y (t) = - 4 (63) 

z (*) = 7772 (A* " J 4 + «i* + (64) 
36 (/j) V 7 

The Hamiltonian constrain gives E = Z\I\ where 2o,i are constants and the singularity condition 
results in the constrain 

_i f (4) 1 + 20 = 0. (65, 

6.2. Power law model W with \i = | 

In this case the Lagrangian eq. ()10p is written as 

21a o 21 a 2 ■ 1 , 7 

L = r d 2 - — — di2- -a 3 i?s. 66 

4i?s 16 Rs 8 

Changing now the variables from (a, i?) to (p, <r) via the relations: 

' 12a 



21\ 3 , „ e 



pe» fl-=— . (67) 



The Lagrangian (|94|) and the Hamiltonian (jll|) become 

t 1 -2 2*2 i t /■ 

L= 2 P "2^ +F °7 2 - 

£ = ^-^* 2 -^. (69) 
where Vo = — A. The Euler-Lagrange equations provide the following equations of motion: 

12ct 

p + /9a 2 + 2Vo^r = (70) 
/-* 

2 e 12CT 
a + -dp+12y ^ = 0. (71) 



The Noether integrals (|43|) , (|44p and the Ermakov-Lewis invariant H5] in the coordinate system 
{u, v} are 



J 5 = 2tE - pp 



1 



I 6 = t 2 £-^p+-p 2 . 



(72) 
(73) 



Z = p 4 & 2 + 4V e 12a . (74) 
Using the Ermakov-Lewis Invariant, the Hamiltonian (|68|) and equation (|7U|) are written: 



1 2 1 S 

-p 2 7 = -E 

2 F 2p 2 

S 

p+ — = 0. 



And the analytical solution of the system is 



p{t)= P2t 2 + Pit + 



(pi) -4E 

4/92 



(75) 
(76) 



(77) 









^tanh 







(tanh o - nP2V / S — 6 arctan h ( ^ 2 't-^ 1 I I — 1 

V L V 2^s yjy 



12 



(78) 



where B (t) = \ h P2 ^ Pl J an d pi,2 > Co are constants with Hamiltonian constrain E = \pi ■ The 

singularity constrain gives (pi) 2 = 4S 

In the case E = the analytical solution is 



p(t)= [ p 2 t 2 + Pit + 



HPi) 



2 \ 2 



exp a (t) 



2 P2 J 
(2p 2 t + Pl ) 



2A^Vo {4,a p 2 2 t + 2a p 2 pi - l) 
The singularity constrain gives pi = 0, then the solution is 



(79) 



(80) 



a(t) 



aott 



(a 2 t- 1)6 



(81) 



In contrast with the claim of [47J this model is analytically solvable and there exists models 
which admit Noether integrals with time dependent gauge functions. 

6.3. A bc CDM model with (6,c) = (1, §) 

Inserting f(R) = (R- 2A) 3 / 2 into eq. (|10p we obtain 

9a 2 - 3 



L = 9aVR - 2Ad 2 + 



+y(^ + 4A)v^r2A 



(82) 



Changing now the variables from (a, R) to (x, y) via the relations: 

9\~^ v 2 

y/£ R = 2A + — (83) 
2/ x 

the Lagrangian ()82p and the Hamiltonian (jlip become 

L = xy + Vq (y 3 + mxy) (84) 

E = xy — Vq (y 3 + fhxy) (85) 

where Vq = g and m = 6A. 

The equations of motion, using the Euler-Lagrange equations, in the new coordinate system are 

x - 3Vb2/ 2 - mV x = (86) 

y - mV y = 0. (87) 
The Noether integrals (|47p . (|48p in the coordinate system {x, y} are 

F 1 = f ?*y-u> ( ?* y (88) 
F^e-^y + ue-^y. (89) 



where oj = y^2A/3. From these we construct the time independent first integral 

<S> = hI 2 = y 2 -0J 2 y 2 . (90) 

The constants of integration are further constrained by the condition that at the singularity 
(t = 0), the scale factor has to be exactly zero, that is, x(0) = 0. 
The general solution of the system (|86[) -(|8T |) is: 

y (t) = A e "* - —e^ (91) 

x (t) = x 1G e ut + x 2G e~" t + - (J 2 e wi + / ie - wi ) 2 + -^L (92) 

The Hamiltonian constrain gives E = oj {x\gI\ ~ X 2G^2) where x\g,2G are constants and the 
singularity condition results in the constrain 

XlG + X 2G + (h + hf + = °- ( 93 ) 

4mw z mu) z 

At late enough times the solution becomes a 2 (t) oc e 2ut 
6.4. Ab c CDM model with (b, c) = (1, f ) 

In this case the Lagrangian eq. ljlOp of the f(R) = (R — 2A) 7 / 8 model is written as 

21a -2 21 ^ ■ £, 1 3 (i2-16A) 
L = -a z w aR a - . (94) 

4(i?-2A)s 16 (i?-2A)s 8 (R-2KY 
Changing now the variables from (a, R) to (u, v) via the relations: 



)« 



The Lagrangian (|94p and the Hamiltonian (jlip become 



L = —u 
2 



1 -2 



m 



A2v 



u z v z + V -u z + V - 



7t z 



£ = -ir — -u"V - Vo—vr - V — 
2 2 4 u 2 

where m = — 28A , Vq = — ^. 

The Euler-Lagrange equations provide the following equations of motion: 



(96) 
(97) 



2 Vom e 
u + uv — u + 2Vq 



I2v 



2 e 
v + -uv + YIVq- 



12v 







0. 



(98) 
(99) 



The Noether integrals (|5i p . (|52p and the Ermakov-Lewis invariant (|53p in the coordinate system 
{u, v} are 



- e 2M E _ e 2M u{i + Ag 2At n 2 

A 

I e -2A tjB _ e -2A tu{( + Ae -2A V _ 
A 

«/, = n V + 4Vbe 12t '. 



(100) 
(101) 

(102) 



where A = §A. 



Using the Ermakov-Lewis Invariant (|102p . the Hamiltonian ()97p and equation (j98|) are written: 



u 2 - Vn^u* 



2« 2 



VqIJI 

u — u + 



E 



0. 



The solution of (|104p has been given by Pinney [19] and it is the following: 

u (t) = ( Ul e 2Xt + u 2 e- 2Xt + 2u 3 



(103) 
(104) 

(105) 



where iti_3- From the Hamiltonian constrain (|103p and the Noether Integrals (|100p . (|10ip we 
find 

E = -2Xu 3 , I + = 2\u 2 , I- = 2\ui. 
Replacing ()105p in the Ermakov-Lewis Invariant (|102p and assuming i^/Owe find: 



exp 



( v ( t )) = 2 i^e- A W (4U + e~ 12A W) S 



where 



^4 (i) = arctan —= (u\e 2Xt + */■; 



+ 4A 2 ui V / 0- 



Then the solution is 

a 2 (t) = 24^ e -^) (V + e- 12A W)" 



^n ie 2A * + U2 e- 2A * + 2n 3 



(106) 
(107) 

(108) 



where from the singularity condition x (0) = we have the constrain u\ + U2 + 2«3 = , or 

2E - (J + + 1_) = 0. (109) 

At late enough time we find A (t) ~ Aq, which implies a?(t) oc e A< . 

In the case where = equations (|1Q3[) . ()1Q4|) describe the hyperbolic oscillator and the 
solution is 

u(t) = sinhAt , 2E = A 2 . (110) 
Prom the Ermakov-Lewis Invariant we have 

.... / AsinhAt \ 6 

y Afi smh At - 12y / |Vb|e" iAt / 

where ui is a constant. The analytical solution is 



9/x / Asinh 7 At \ 6 

a (*) = == (112 

^Awi sinh At - \2yJ\V^\e- 2Xt ) 

7. Noether symmetries in spatially non flat f(R) models 

In this section we study further the Noether symmetries in non flat f(R) cosmological models. 
In the context of a FRW spacetime the Lagrangian of the overall dynamical problem and the 
Ricci scalar are 

L = Qfaa 2 + 6f"Ra 2 a + a 3 (f'R - f) - 6Kaf (113) 

fii a 2 + K\ 
R = 6(-+ \ 114 
\a a z J 

where K is the spatial curvature. Note that the two dimensional metric is given by eq. (|35p while 
the "potential" in the Lagrangian takes the form 

V K (a,R) = -a 3 (f'R-f) + Kaf. (115) 

Based on the above equations and using the theoretical formulation presented in section [SJ we find 
that the f(R) models which admit non trivial Noether symmetries are the f(R) = (R — 2A) 3 / 2 , 
f(R) = R 3 / 2 and f(R) = R 2 . The Noether symmetries can be found in section 

In particular, inserting f(R) = (R — 2A) 3 / 2 into the Lagrangian (|113p and changing the 
variables from (a,R) to (x, y) [see section IS\3"] we find 

L = xy + Vq (y 3 + fhxy) — Ky (116) 

E = xy - V (y 3 + mxy) + Ky (117) 
where K = 3(6 1 / 3 AT). Therefore, the equations of motion are 

x - 3V y 2 - fhV x + K = 
y - mV y = . 



The constant term K appearing into the first equation of motion is not expected to affect the 
Noether symmetries (or the integrals of motion). Indeed we find that the corresponding Noether 



symmetries coincide with those of the spatially flat f(R) = (R — 2A) 3 / 2 model. However, in the 
case of K ^ (or the analytical solution for the x-variable is written as 

K 

x K {t) = x(t) + (118) 

where x(t) is the solution of the flat model K = (see section loT3|) . Note that the solution of 
the y-variable remains unaltered. 

Similarly, for the case of the f(R) = R 3 ^ 2 model the analytical solution is 

z K (t) = z{t) + K (119) 

where z (t) is the solution of the spatially flat model (see section UTTj) . 



8. Conclusion 

In the literature the functional forms of f(R) of the modified f(R) gravity models are mainly 
defined on a phenomenological basis. In this article we use the Noether symmetry approach to 
constrain these models with the aim to utilize the existence of non-trivial Noether symmetries 
as a selection criterion that can distinguish the f(R) models on a more fundamental level. 
Furthermore the resulting Noether integrals can be used to provide analytic solutions. 

In the context of f(R) models, the system of the modified field equations is equivalent to 
a two dimensional dynamical system moving in M 2 (mini superspace) under the constraint 
E =constant. Following the general methodology of [31,38], we require that the two dimensional 
system admits extra Noether symmetries. This requirement fixes the / (R) function and the 
analytical solutions are computed. It is interesting that two well known dynamical systems 
appear: the anharmonic oscillator and the Ermakov-Pinney system. We recall that the field 
equations of the A— cosmology is equivalent with that of the hyperbolic oscillator. 
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Appendix A. Special solutions for the Power law model R n 

The case / (R) = R n is not Liouville integrable via Noether point symmetries. The 
zero order invariant will be used in order to find special solutions. Inserting f(R) = 
R n , (n / 0, 1, §, |) into eq.([10]) we obtain 

L (a, a, R, Rj = QnaR n ~ l a 2 + 6n (n - 1) a 2 R n ~ 2 aR + (n - 1) a 3 R n (A.l) 
and the modified field equations are 

a + -d 2 --aR = (A.2) 
a 6 

R + ^R 2 - -^-r*d 2 + -aR - } n ~\ R 2 = (A.3) 
R n-lr a bn[n — \) 

E = GnaR^a 2 + 6n (n - 1) a 2 R n ~ 2 aR - (n - 1) a 3 R n . (A.4) 



The Noether symmetry (|54p is also and a Lie symmetry, hence we have the zero order 
invariants 

a = at~ N , Rq = Rt~ 2 . (A.5) 



Applying the zero order invariants in the field equations (|A.2p - (|A.4p and in the Noether 
integral (|55p we have the following results. 

The dynamical system admit a special solution of the form 



a = a t N , R = 6N(2N -l)t~ 2 (A.6) 
where the constants N, E and J7 are 

N = - , E = , J 7 = 



2 



or 



or 

2 

= -77. 

3 



(2n - 1) (n - 1) 

A r = ^ , E = , J 7 = 

n — z 

£ = (4n - 3)™" 1 (I3n - 8n 2 - 3) 4 , I 7 = 0. 



Another special solution is the deSitter solution for n = 2 

a = a e Hot , R = 12H$ (A.7) 
where I-j = and the spacetime is empty i.e. E = 0. 
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